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We  have  studied  grid-generated  turbulence  in  a  water  tunnel  at  moderate  Reynolds 
numbers.  The  method  used  was  the  standard  one  of  Laser  Doppler  Velocimetry  (LDV) 
and  an  novel  scheme,  which  we  call  photon  homodyne  correlation  spectroscopy(HCS). 
With  LDV,  we  measured  the  probability  density  function  of  velocity  differences  P{Sv{()) 
on  varying  spatial  scales  £,  by  invoking  the  frozen  turbulence  hypothesis.  The  HCS 
technique  permits  measuring  P  without  using  this  hypothesis.  Of  special  interest  to 
us  was  the  behavior  of  the  system  at  and  above  a  Reynolds  number  (Rcc)  where  the 
turbulence  becomes  self-similar,  in  that  <  6v{i)^  >~  £^.  Above  Rcc  the  exponent  C 
increases  from  0  to  2/3  with  increasing  Re, 
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This  research  was  a  study  of  grid-generated  turbulence  by  a  novel  method  of  photon  homodyne 
correlation  spectroscopy  (HCS).  The  scheme  enables  one  to  measure  velocity  differences  V{R)  on  a 
spatial  scale  R,  with  no  need  to  invoke  the  Taylor’s  frozen  turbulence  hypothesis.  In  this  work  we 
compared  such  measurements  with  with  Laser  Doppler  Velocimetry  measurements  of  the  velocity 
at  a  point,  v{t).  From  the  latter  experiment,  we  deduced  the  velocity  difference  V{R)  with  the 
Taylor  assumption,  R  =  6t/U,  where  U  is  the  mean  flow  velocity  and  6t  is  the  time  interval  between 
velocity  measurements. 

In  that  work  it  was  observed  that  the  time-time  averaged  velocity  difference,  <  V{R)  >,  was  of 
self-similar  form,  <  >=  aR^,  where  C  is  a  function  of  Re,  and  o  is  a  constant,  which  was  also 

measured.  The  AFOSR  support  enabled  us  to  build  a  much  larger  tunnel  and  to  span  a  wider  range 
of  Reynolds  numbers.  Before  summarizing  our  new  results,  I  describe  the  HCS  measuring  scheme 
in  elementary  terms. 

If  a  light  beam  strikes  a  single  particle  moving  with  velocity  v,  the  frequency  of  the  light  is  shifted  by 
an  amount  k  -  v,  where  k  is  the  the  scattering  vector  of  the  light.  If  the  scattering  angle  with  respect 
to  the  incident  beam  is  6,  then  the  magnitude  of  k  is  {4ir/X)ain(9/2),  where  A  is  the  wavelength  of 
light.  This  frequency  shift  is  the  basis  for  laser  Doppler  Velocimetry  (LDV).  If  now,  the  incident 
beam  strikes  two  particles,  moving  with  velocities  vt  and  vj,  and  if  the  scattered  light  falls  on  a 
photomultiplier,  the  output  of  this  device  will  be  a  photo-current  modulated  at  the  beat  frequency 
k  vi  -  kvj).  Thus,  one  has  a  measurement  of  the  difference  in  velocity  of  the  two  particles.  This 
is  the  basis  of  HCS. 

It  turns  out  that  one  can  measure  velocity  differences  even  when  the  beam  traversing  the  turbulent 
fluid  is  scattered  by  a  large  number  of  particles.  In  this  case,  the  beatings  of  Doppler  shifts  from 
many  particle  pairs,  causes  the  scattered  intensity,  I(t),  to  fluctuate  in  time.  These  temporal  fluctu¬ 
ations  can  be  measured  by  a  frequency  analyzer,  or  more  conveniently,  by  measuring  the  intensity 
correlation  function,  <  >,  where  the  brackets  designate  an  average  over  f'.  If  one  images 

the  scattered  light  beam  on  a  slit  of  adjustable  width  L,  then  the  (normalized)  intensity  correlation 
function  <  >  /  <  /(<')  >^>=  5(f).  has  the  form. 


ff(t)  =  I  -b  /  i(fl)d(«)  /  P(V(R))cos(kV(R))dV(R)  =  1  -h  C(t), 


where  6(i2)  =  (2/£)(l  —  RfL)  is  the  number  fraction  of  seed  particles  separated  by  a  distance  R, 
in  the  slit,  and  P{R)  is  the  probability  density  function  for  velocity  differences  on  a  scale  {R). 
The  symbol,  V{R)  denotes  the  projection  of  the  velocity  difference  V(/Z)  along  the  direction  of  the 
scattering  vector  k. 

We  note  severad  features  of  the  above  equation.  The  cosine  factor  is  just  the  beat  frequency  associated 
with  velocity  differences  V’(/2),  coming  from  particle  piurs.  Because  the  speed  of  light  is  so  very  large, 
the  experiment  effectively  measures  simultaneously,  velocity  difference  from  the  particles  constituting 
a  pair  (It  turns  out  that  the  scattering  from  triplets  of  particles,  etc,  average  to  zero).  Hence  the 
Taylor  hypothesis  need  not  be  invoked.  In  principle  a  knowledge  of  P(R)  is  sufficient  to  determine 
all  the  moments  of  V(R),  from  which  the  spectrum  of  fractal  exponents  can  be  extracted  -  if  indeed 
the  turbulent  energy  lies  on  a  fractal  or  multifractal.  In  practice,  one  cannot  measure  the  higher 
moments  of  P(R)  with  the  HCS  scheme,  but  one  can  learn  a  lot  about  the  functional  form  of  P(R). 
Because  the  direction  of  the  incident  and  scattered  beams  eire  at  the  experimenters  disposal,  velocity 
differences  V(/Z)  can  be  measured  in  any  direction,  so  the  presence  or  absence  of  isotropy  of  the 
turbulence  can  be  determined.  Like  LDV,  the  HCS  scheme  is  non-invasive. 

This  grant  enabled  us  to  construct  a  rather  fine  tunnel,  roughly  10  cm  x  10  cm  in  cross  section.  An 
important  feature  of  this  tunnel  was  that  its  flat  walls  were  of  optical  quality,  so  that  the  scattering 
vector  was  well  defined,  and  a  sharp  image  could  be  formed  on  the  slit.  The  mesh  size  M  of  the  grid 
in  this  new  tunnel  was  0.85  cm.  The  grant  also  enabled  me  to  buy  a  laser  Doppler  velocimeter,  so 
that  it  was  possible  to  compare  the  results  of  the  HCS  and  LDV  noeasurements  at  the  same  point 
in  the  turbulent  stream,  namely  21  cm  down  stream  from  the  grid. 

In  the  LDV  measurements  we  concentrated  on  measuring  the  probability  density  of  velocity  differ¬ 
ences  P‘(V{R)),  so  that  this  function  could  be  compared  with  the  HCS  measurements  of  P(V(R)) 
introduced  above.  If  the  Taylor  hypothesis  were  applicable,  P  and  P'  would  be  equal,  assuming,  of 
course  that  the  values  of  R  were  the  same  in  both  experiments.  It  turns  out  that  the  HCS  scheme 
is  most  suited  for  studying  velocity  differences  in  small  scales,  R,  whereas  the  LDV  scheme  yielded 
the  most  accurate  information  at  large  values  of  R.  Nevertheless,  the  overlap  interval  of  the  two 
measuring  schemes  was  significant,  namely  0.5  mm  <  R  <  2.5  mm.  With  the  present  optics  we 
probe  eddies  as  small  as  0.1  mm  with  HCS,  the  limitation  being  the  width  to  which  we  can  focus 
the  beam  segment  which  is  imaged  on  the  slit. 

We  have  found  that  if  Re  exceeds  some  critical  value,  which  we  call  Rce,  the  probability  density 
P(V(R))  has  the  simple  scaling  form,  P(V(R)/«(R)),  where  u(R)  is  the  characteristic  width  of  P. 
Over  a  decade  in  R,  this  width  varies  algebraically  with  R,  i.e.,  u(R)  ~  In  the  Kolmogorov 

theory,  the  exponent  (  =  1/3.  We  find,  on  the  other  hand,  that  this  exponent  increases  from 
approximately  zero  to  1/3  as  Re  is  increased  above  Rcg.  At  our  largest  achievable  values  of  Re,  C  — 
1/3.  One  interpretation  of  this  observation  is  that  the  turbulence  is  increasing  from  sheet-like  (fractal 
dimension  D  =  2)  to  volume-filling  (D=3)  as  Re  is  increased.  In  quite  independent  experiments, 
carried  out  in  a  different  way,  Sreenivasan  and  his  students  have  come  to  the  same  conclusion. 

Prom  a  measurement  of  G{t)  one  can,  in  principle,  deduce  P(V')  using  Eq.  1.  One  sees  from  that 
equation  that  P  is  closely  related  to  the  Fourier  transform  of  G(().  In  practice  it  is  often  difficult 
to  determine  a  function  from  its  experimentally  measured  Fourier  transform,  and  thus  it  was  here. 
All  we  can  say  is  that  we  can  recover  the  measured  G(()  with  a  P  that  is  the  product  of  Gaussian 
and  Lorentzian  factors,  i.e. 

P(R)  =  const  e*p(-V(r)/2u'(R)*)  x  1/(V(R)’  +  u(R)’),  (2) 

with  u'(R)  roughly  five  times  u(R).  This  large  ratio  implies  that,  in  the  range  of  our  measurements, 
P  departs  only  slightly  from  Lorentzian  form.  The  LDV  measurements  are  also  consistent  with  a 
P(R)  of  the  above  form,  but  in  this  case  u  and  u'  are  almost  equal.  The  measurements  leave  no  doubt 
that,  in  the  domain  of  eddies  sizes  R  and  Reynolds  numbers  probed  by  these  experiments,  the  time 
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domain  measurements  (LDV)  and  the  space  domain  measurements  (HCS)  yield  very  different  forms 
for  P{V).  Almost  all  existing  measurements  of  this  probability  density  are  made  in  the  time  domain, 
with  the  spatial  information  extracted  by  invoking  the  Taylor  hypothesis.  Recently  Sreenivasan  has 
taken  some  measurements,  using  a  single  hot  wire  annemometer  and  also  a  pair  of  them.  He  too 
finds  a  failure  of  the  frozen  turbulence  assumption.  This  past  week  I  will  give  an  invited  colloquium 
at  Yade  and  had  an  opportunity  to  compare  notes  with  him. 

I  have  summarized  the  work  completed  under  the  AFOSR  grant.  The  experiments  in  our  water 
t  tin  net  are  continuing,  and  we  are  also  probing  the  transition  from  chaos  to  turbulence  in  a  couette 
cell.  We  are  also  exploring  the  possibility  of  measuring  the  probability  density  for  the  vorticity, 
P(u(R))  by  optical  means.  The  completed  work  has  resulted  in  a  number  of  publications,  as  you  see 
in  the  addendum.  In  addition  I  have  been  invited  to  present  our  results  at  various  meetings,  both 
in  the  US  and  abroad.  Professor  Sirivat,  Prfrfessor  Tong,  Mr.  Pak  and  I  are  very  grateful  for  the 
support  we  have  received  from  the  AFOSR.  The  work  done  under  this  grant  constitutes  the  PhD 
thesis  of  Mr.  Pak. 
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List  of  of  publications  that  came  from  the  AFOSR  support: 

1.  “A  Light  Scattering  Study  of  Turbulence”,  (W.  1.  Goldburg,  P.  Tong,  and  H.  K.  Pak), 
Physica  D  3S  (1989)  134-140. 

2.  “Scaling  Laws  in  Weak  Turbulence”,  (with  H.  K.  Pak,  W.  1.  Goldburg  and  P.  Tong), 
(in  press). 

3.  “An  Experimental  Study  of  Weak  Turbulence”,  (with  H.  K.  Pak  and  W.  1.  Goldburg), 
“Fluid  Dynamics  Research”  (in  press). 

4.  "Measuring  the  Probability  Distribution  of  the  Relative  Velocities  in  Grid-Generated 
Turbulence”  H.K.  Pak,  and  W.I.  Goldburg,  (submitted  for  publication). 

Invited  talks  that  are  an  outgrowth  of  this  work: 

1.  International  Conference  on  Fractals  in  Physics 

(A  conference  in  honor  of  B.  Mandelbrot),  Nice,  1989. 

2.  SIAM  Conference  on  Turbulence,  Orlando,  May,  1990. 

3.  Gordon  Conference  on  Fractals,  Plymouth  NH,  August,  1989. 

4.  International  Workshop  on  Novel  Experiments  and  Data  Processing 
for  Basic  Understanding  of  Turbulence,  Tokyo,  October,  1990. 

5.  Program  on  Spatially  Extended  Nonequilibrium  Systems,  Santa  Barabara, 

Aug.  1992. 
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A  Light  Scattering  Study  of  Turbulence 
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A  LIGHT  SCATTERING  STUDY  OF  TURBULENCE 
W.I.  GOLDBURG,  P.  TONG  '  and  H.K.  PAR 

Department  ni  Physics  and  Aslronomv.  i'niversity  of  Pittsburgh.  Pittsburgh.  P.-i  15260.  L'S.4 


By  scattering  light  from  a  turbulent  tluid  seeded  with  small  particles,  one  obtains  information  about  turbulent  velocity  fluctua¬ 
tions  over  varying  spatial  scales.  R.  The  measured  intensity  autocorrelation  function,  git),  is  related  to  the  probability  density 
P( !  ( /? ) )  of  finding  velocity  fluctuations  of  magnitude  f  '(  R )  associated  with  eddies  of  size  R.  The  measurements  described  here 
strongly  suggest  that  the  energy<ontaining  eddies  occupy  a  fractal  region  whose  dimension  (or  spectrum  of  dimenstons)  increases 
with  the  Reynolds  number  Re  when  Re  exceeds  some  threshold  value. 


I.  Introduction 

In  his  seminal  book.  The  Fractal  Geometry  of  Na¬ 
ture.  Benoit  Mandelbrot  ( 1  ]  makes  clear  his  deep  in¬ 
terest  in  the  geometrical  nature  of  turbulence.  As  he 
points  out.  the  description  of  the  visual  appearance 
of  a  turbulent  fluid,  such  as  smoke  curling  up  from  a 
cigarette,  taxes  our  powers  of  description.  It  seems 
that  present-day  speech  is  not  well  suited  to  evoking 
the  image  of  self-similar  structures.  After  all.  it  takes 
a  series  of  images,  one  magnified  with  respect  to  the 
other,  to  identify  fractal  structures.  And  turbulence 
IS.  by  all  evidence,  a  fractal  thing  at  its  roots[2]. 

There  are  many  ways  of  revealing  the  fractal  or 
spotty  nature  of  a  turbulent  fluid.  One  technique  is  to 
measure  the  time  variation  of  the  square  of  the  veloc¬ 
ity  at  a  point  in  the  fluid  [3].  Another  is  to  add  a 
small  amount  of  long-chain  molecules  to  the  fluid  and 
observe  it  through  crossed  polaroids  [4].  The  mole¬ 
cules  are  locally  aligned  by  turbulent  shear  forces. 
These  molecules,  being  anisotropic  scatterers.  depo¬ 
larize  the  light  in  regions  where  the  shear  is  large, 
making  the  local  structure  of  the  strong  vorticity  di¬ 
rectly  visible. 

Herein  we  describe  expenments.  carried  out  at  the 
University  of  Pittsburgh,  which  provide  a  new  ap- 
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proach  to  the  study  of  the  small-scale  structure  of  tur¬ 
bulence.  The  method  involves  a  measurement  of  the 
autocorrelation  function  of  the  light  intensity  scat¬ 
tered  by  small  particles  suspended  in  the  turbulent 
fluid.  For  this  technique  there  is  no  need  to  invoke 
the  “frozen  turbulence  assumption"  to  translate  tem¬ 
poral  information  to  spatial  information.  According 
to  this  assumption,  small-scale  eddies  (the  ones  of 
interest),  are  transported  past  a  velocity  measuring 
device  with  the  mean  velocity  C^of  the  flow.  If  these 
small-scale  eddies  remain  intact  for  a  long  enough 
time,  a  time  record  of  the  velocity  v(  f )  at  a  point  will 
reveal  spatial  features  of  the  flow  through  the  equa¬ 
tion  v{t)  =  v{x/U).  The  frozen  turbulence  assump¬ 
tion  fails  unless  the  velocity  fluctuations  i'(R)  asso¬ 
ciated  with  eddies  of  size  R  are  uncoupled  from  the 
larger-scale  eddies. 

The  technique  of  photon  correlation  homodyne 
spectroscopy  (HS)  [5].  which  we  have  used  in  our 
experiments,  is  that  of  recording  the  beating  of  scat¬ 
tered  light  waves  that  have  been  Doppler  shifted  by 
pairs  of  particles  seeded  in  the  turbulent  fluid.  The 
technique  was  introduced  many  years  ago  by  Bourke 
et  al.  [6],  but  seems  largely  to  have  been  ignored. 
Being  an  optical  technique,  it  permits  non-invasive 
observation  of  velocity  fluctuations  at  very  small 
scales. 

The  homodyne  scheme  is  readily  understood  from 
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fig.  1,  which  shows  two  moving  particles  at  a  partic¬ 
ular  instant  of  time  when  their  separation  is  R  and 
their  velocities  are  r,  and  r,.  The  seed  particles  are 
small  enough  that  they  scatter  light  isotropically.  A 
photodetector  ( PMT),  located  at  an  angle  6  with  re¬ 
spect  to  the  incident  beam,  receives  the  light  from 
both  particles.  The  scattered  light  from  each  particle 
is  Doppler  shifted  by  an  amount  At,  and  respec¬ 

tively,  where  A  is  the  scattering  vector,  of  magnitude 
A=(4nn/A)  sin(J^).  Here  X  is  the  vacuum  wave¬ 
length  of  the  light  (A  =  488  nm  in  our  experiments), 
and  n  is  the  refractive  index  of  the  turbulent  fluid:  in 
our  case  the  fluid  was  water.  The  photomultiplier, 
which  receives  the  light  from  the  particle  pair,  is  a 
square-law  detector,  so  that  its  output  current,  /(r), 
contains  a  beating  term  proportional  to 
cos(AKj(/?)r],  where  V/,  is  the  projection  of  the  ve¬ 
locity  difference  Vi-Vz  along  the  direction  of  A. 
Henceforth  the  subscript  on  will  be  dropped, 
but  its  R  dependence  will  be  retained. 

The  essential  aspect  of  turbulence  is  that  the  veloc¬ 
ity  difference  between  two  points  in  the  fluid  de¬ 
pends  on  the  separation  R  of  these  two  points.  Ac¬ 
cording  to  the  theory  of  Kolmogorov  [6),  the 
moments  of  the  velocity  fluctuations  obey  a 
scaling  law 


iViRyy^uiRy-R^ .  (1) 

with  Jn.  The  homodyne  technique  is  well  suited 
to  measure  the  lower  moments  of  V^R),  but  not  the 
higher  moments.  On  the  other  hand,  the  method 
yields  information  about  the  functional  form  of  the 
probability  density  P(  K(/{) ),  that  two  points  in  the 
fluid,  separated  by  a  distance  R,  have  velocity  differ¬ 
ence  lying  within  K(R)  and  K(/{)-t-dK(/i). Our  cen¬ 
tral  finding  is  that  P(y{R))  is  well  represented  by  a 
Lorentzian  function, 

P(V(R))^{\Ar[V(R)/u(R)]-}-' ,  (1) 

for  relatively  small  values  of  K(/? ).  We  also  find  that 
the  scaling  velocity  u{R)~R^,  where  C  is  a  function 
of  Reynolds  number.  Re.  The  measurements  were 
made  at  very  modest  values  of  the  Reynolds  number. 
In  fact  the  turbulence  was  so  weak  that  one  might  not 
have  thought  the  flow  would  exhibit  the  self-similar¬ 
ity  which  was  indeed  observed.  Throughout  this  pa- 
-  per  the  Reynolds  number  is  defined  as  Re=Uto/v. 
where  U  is  the  mean  velocity  of  the  flow,  /q  is  the  outer 
scale  of  the  turbulence,  and  v  is  the  kinematic  viscos¬ 
ity  of  the  fluid. 


Fig.  I  A  schematic  diagram  showing  scattering  geometry.  The  scattering  vector  *=*,-*o.  where  *,  and  ko  are  the  scattered  and  incident 
wave  vectors  respectively. 
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2.  Experimental 

The  detailed  experimental  setup  can  be  found  in 
ref.  [  8  ] .  The  fluid  flow  was  generated  in  a  closed 
water  tunnel  comprised  of  a  cylindrical  pipe  and  a 
pump  of  variable  speed.  The  turbulence  is  generated 
by  a  grid  within  the  pipe.  The  grid  can  be  removed  to 
permit  study  of  wall-generated  turbulence  (pipe 
flow ) .  A  baffle  section  placed  in  the  high-pressure  side 
of  the  grid,  suppresses  the  turbulence  generated  by 
the  pump  and  by  those  sections  of  pipe  on  the  high- 
pressure  side  of  the  grid.  With  this  arrangement  all  of 
the  turbulence  is  generated  by  the  grid  only.  In  most 
of  the  experiments  discussed  here  the  diameter  of  the 
pipe  was  4.4  cm.  and  the  aperture  size  of  the  grid  was 
3. 1  mm.  These  parameters  are  taken  to  be  /(>  in  cal¬ 
culating  the  Reynolds  number.  The  measurements 
were  made  28  cm  downstream  from  the  grid.  The 
water  which  flowed  through  the  pipe  was  seeded  with 
polystyrene  spheres  60  nm  in  diameter.  These  parti¬ 
cles  were  small  enough  to  scatter  light  isotropically 
and  tn  sufficient  concentration  that  their  mean  sepa¬ 
ration  was  much  less  than  the  Kolmogorov  dissipa¬ 
tion  length  4|.  which  was  estimated  to  be  a  fraction  of 
a  millimeter. 

On  the  downstream  side  of  the  grid  there  is  an  op¬ 
tically  transparent  section  of  piping  to  admit  the  in¬ 
cident  laser  beam  and  observe  the  scattering.  Be¬ 
cause  the  flow  is  seeded,  a  thin  column  of  the  scattered 
light  is  produced  in  the  water  and  that  light  is  imaged 
with  a  lens,  on  a  slit  of  vanable  width.  L.  By  varying 
L.  the  homodyne  scheme  permits  the  probing  of  ve¬ 
locity  fluctuations  !'(/?)  from  the  smallest  scale  /j  to 
that  of  the  width  of  the  slit.  L. 

Using  a  standard  light  scattering  apparatus  and  a 
digital  correlator,  we  measure  the  intensity  autocor¬ 
relation  function.  ^(/)  =  </(f' )  1(1'  )>^ 

where  l(t)  is  the  scattered  light  intensity  measured 
at  scattering  angle  9.  and  the  angle  brackets  represent 
a  time  average  over  t  .  One  can  show  (8)  that  the 
correlation  function  g(  i )  has  the  following  form: 

?(f)  =  l  (7(0  .  (3) 

The  geometncal  factor  /  (.■( )  is  of  order  unity  if  the 


photodetector  receives  light  from  only  one  coherence 
area  [5].  All  of  the  interesting  physics  is  contained 
in  G(t).  which  is  proportional  to  a  sum  of  the  time- 
averaged  phase  factors  cos(^7T)  coming  from  the 
Doppler  shift  of  all  particle  pairs  in  the  scattering 
volume.  The  function  G(i)  can  be  written  as  [8] 

L  X 

G(0=  jd/?h(R)  j  dy(R)  P(y(/i}) 

0  -  X 

.  Xcos[/cT(R)0  .  (4) 

where  7i(/?)  is  the  probability  of  finding  a  particle 
pair,  separated  by  R.  in  the  columnar  region  of  length 
L. 

If  the  image  on  the  slit  is  taken  to  be  quasi-one- 
dimensional.  which  is  valid  when  the  slit  width  re¬ 
mains  large  compared  to  the  diameter  of  the  laser 
beam.  h(R)=i2{  \~R/L)/L.  Note  that  the  inner  in¬ 
tegral  in  eq.  (4)  is  the  Fourier  cosine  transform  of 
P(f''(R)}.  and  the  (7  ( r )  may  be  thought  of  as  a  trans¬ 
form  of  the  characteristic  function.  If  the  probabil¬ 
ity  density  P(U(R))  has  the  scaling  form 
P(  y(R))=Q[V(R)/u(R)]/u(R),cq.  (4)  becomes 

L 

G(i)=  ^  dRhiR)  F{kiu{R))  .  (5) 

where  F  is  the  Fourier  cosine  transform  of  (?[  U(/?)/ 
u(R)].  It  is  easy  to  show  that  the  scaling  law  in  eq. 

( 1 )  follows  if  the  probability  density  function 
F(  U(  /? ) )  has  the  above  mentioned  scaling  form. 

The  above  equations  for  g(l)  have  quite  general 
validity.  They  hold,  for  example,  even  if  the  fluid  is 
stationary,  and  the  seed  particles  are  undergoing 
Brownian  motion  only.  In  that  case.  I'(R)  is  inde¬ 
pendent  of  R  and  P( !'(/?))  is  a  Gaussian  function. 
Then  the  function  G(t)  is  an  exponentially  decaying 
function  [5],  G(t)=exp(-2Dk-i).  where  D  is  the 
diffusivity  of  the  Brownian  particles,  and  is  given  by 
Stokes'  law.  This  contribution  to  the  decay  of  Git) 
will  be  present,  even  when  the  fluid  is  turbulent. 
However,  in  a  turbulent  fluid.  decay  time  T  of 
G{t)  is  much  shorter  than  the  diffusive  decay  time. 
Ta=  \/2Dk-,  so  that  the  latter  contribution  can  be 
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safely  ignored.  From  eq.  (  5 ).  it  follows  that  the  tur¬ 
bulent  decay  time  should  be  of  the  order  of  /■= 
\  lku{L).  because  the  fastest  decay  rate  is  associated 
w  ith  the  largest  eddies  of  size  L. 


3.  Results 

Over  a  wide  range  of  slit  widths  and  Reynolds 
numbers,  we  find  that  the  function  (/(/ )  exhibits  the 
scaling  form. 

(>(  kt.  L.  Re )  =  f/ ( V  (  .  ( 6 ) 

with  K^k  'L'i.  This  scaling  behavior  of  0'(t)  is  ob¬ 
served  only  when  the  Reynolds  number  exceeded  a 
certain  value  Re^.  In  the  case  of  the  grid  flow  de¬ 
scribed  above.  Re^  was  roughly  500.  which  corre¬ 
sponds  to  much  weaker  turbulence  than  that  one  nor¬ 
mally  associates  with  scaling  behavior.  I*  is  quite 
possible  that  the  scaling  behavior  is  seen  at  such  small 
values  of  Re  because  the  simultaneous  velocity  dif¬ 
ference  l  (  R )  IS  measured  and  no  frozen  turbulence 
assump..on  is  needed  in  the  data  analysis. 

The  exponent  in  eq.  ( 6 )  was  measured  as  fol¬ 
lows.  For  a  fixed  slit  width  L  and  a  fixed  Re.  (id) 
was  measured  at  several  scattering  angles  and  hence 
several  values  of  k.  .All  of  the  plots  of  log  [(7(f))  vci- 
sus  log(f)  could  be  superimposed  by  horizontal 
translation  of  one  graph  with  respect  to  another.  The 
amount  of  translation.  6{k),  is  found  to  be  roughly 
proponional  to  k,  i.e.  n=\.  when  Re  exceeded  the 
critical  value  RCc.  However,  in  the  absence  of  flow. 
u  =  2.  as  expected  for  Brownian  motion  of  the  seed 
particles.  Similar  measurements  were  made  in  which 
k  and  Re  were  held  fixed,  and  L  was  vaned.  .Again  all 
the  plots  of  log[G(f )  1  versus  log(f )  could  be  super¬ 
imposed.  yielding  the  result  K-^kL't.  as  long  as  Re 
exceeded  Rec.  In  these  experiments.  C  is  found  to  be 
Re  dependent.  We  return  to  this  important  observa¬ 
tion  below. 

Fig.  2.  a  log-log  plot  of  G(k)  versus  k.  shows  the 
scaling  behavior  of  Gd)  discussed  above.  The  mea¬ 
surements  correspond  to  several  values  of  scattering 
angle,  or  /f-value.  several  slit  widths  and  at  vanous 


Fig.  i.  The  scaling  funciion  0'(  k  )  versus  K  =  tfu(  LU  in  pipe  tlow 
and  grid  flow. 

Reynolds  numbers.  In  one  set  of  measurements 
(closed  circles),  the  gnd  was  present;  in  the  other  two 
sets  (crosses  and  triangles),  the  grid  was  removed 
(pipe  flow).  The  correlation  functions  (Sd)  have 
been  horizontally  (and  vertically)  translated  so  that 
they  coincide.  In  the  pipe  flow  measurements,  the 
Reynolds  number  is  based  on  the  pipe  diameter, 
making  it  an  order  of  magnitude  larger  than  that  for 
the  grid  flow,  even  when  the  mean  flow  velocities  L 
arc  comparable  in  both  cases. 

An  alternative  way  to  determine  the  exponent  2  was 
to  plot,  on  a  double  logarithmic  scale,  the  slit-width 
dependence  of  the  decay  time.  T.  of  G(  f ) .  keeping  Re 
and  k  fixed.  As  is  shown  in  fig.  3.  linear  vanation  of 
log(  T)  with  log(L)  was  seen  at  intermediate  values 
of  L.  The  data  in  fig.  3  were  obtained  in  the  grid  flow 
at  three  different  Reynolds  numbers  Re  =  460.  1400. 
and  2200.  Since  7'=sl/A:u(L)  and  u(L)-L-.  the 
slope  of  this  line  yields  the  exponent  2-  which  is  I  /  3 
in  the  Kolmogorov  theory.  We  have  venfied  that  the 
power  law  behavior  at  large  L  was  limited  by  the  outer 
scale,  /n.  of  the  turbulence.  At  small  values  of  L  the 
beam  diameter  was  no  longer  negligibly  small,  w  hich 
could  account  for  the  decrease  in  T at  small  values  of 
L.  Imperfections  in  the  optical  system  may  also  be 
responsible  for  the  decrease  in  T at  small  slit  widths. 

The  behavior  of  r(  L )  at  small  L  has  more  recently 
been  reexamined  in  a  water  tunnel  of  much  superior 
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Fig.  3  The  decay  lime  TlL)  venus  slit  width  L  in  jnd  flow.  The 
number  below  a  line  is  the  slope  of  that  line. 

design  to  that  used  in  the  studies  reported  above.  In 
this  experiment,  the  optically  transparent  pipe,  where 
g(  i )  was  measured,  was  square  in  cross  section,  rather 
than  cylindrical,  so  that  the  laser  beam  was  undis¬ 
torted  in  passing  through  it.  In  this  square  pipe,  the 
beam  diameter  was  less  than  0. 1  mm,  which  is  smaller 
than  the  smallest  value  of  L  at  which  g{t)  was  mea¬ 
sured.  Using  laser  Doppler  velocimetry  and  invoking 
the  frozen  turbulence  assumption  one  can  determine 
the  smallest  eddy  size  U.  At  Res 8 SO,  we  obtained 
'd=0.4  mm.  At  values  of  L  between  0.4  mm  ( =  4) 
and  0. 1  mm.  the  decay  time  of  G(  r )  became  indepen¬ 
dent  of  L,  i.e.  C*0  when  Z,<4.  This  result  is  very 
different  from  the  Kolmogorov  prediction.  C=:  I  when 
L<U. 

From  the  straight-line  segment  (solid  line  in  Tig.  3 } 
we  can  extract  the  slope  C  which  shows  a  Re-depen¬ 
dent  feature.  Fig.  4  shows  C  as  a  funaion  of  Re  for 
both  pipe  flow  and  grid  flow  ( insert ).  The  exponent 
C  is  seen  to  increase  from  0  to  «  I  /3  (the  Kolmogo¬ 
rov  value)  as  Re  is  increased.  When  the  Reynolds 
number  is  below  Re,,  0(t)  fails  to  exhibit  scaling  be¬ 
haviour.  The  measured  Re^ «  3(X)-4(X)  in  the  grid  flow 
and  Rcr  %  3(XX)-4(XX)  when  the  grid  was  removed. 
Measurements  in  the  improved  water  tunnel  give 
similar  results.  These  observations  are  consistent  with 


Fig  4.  The  exponent  C  u  a  function  of  Re  in  pipe  flow.  The  solid 
line  IS  drawn  by  eye  through  the  dau  points,  and  ihe  dashed  curve 
shows  the  oscillatory  behavior  of  v  The  inset  shows  C  versus  Re 
in  gnd  flow. 


Fig  5.  A  plot  of  log[  C(  9(.  Z. )  ]  versus  t  in  pipe  flow  at  indicated 
parameters. 

the  notion  that  the  turbulence  becomes  increasingly 
three-dimensional  as  Re  is  increased  above  RCc  and 
that  in  the  vicinity  of  Rcc,  the  turbulence  is  two-di¬ 
mensional  [9]. 

We  now  turn  to  the  discussion  about  the  functional 
formofC(f).  Fig.  5  is  a  semilog  plot  of  G({)  versusf 
in  pipe  flow  at  the  indicated  values  of  L.  6,  and  Re. 
The  straight  line  is  a  linear  fit  to  the  dau  points  at 
small  {.  It  is  seen  that  only  at  large  time  does  the  curve 
sun  to  deviate  from  the  linear  behavior.  If  we  as¬ 
sume  that  the  charaaeristic  function  F(ku{R)t)  in 
eq.  (5)  has  the  form  f  ~exp(  -fcu(R)/),  G(t)  then 
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becomes  an  incomplete  gamma  function  with  ku(L)t 
as  Its  argument  [8].  This  equation  is  well  fitted  to 
our  measurement  of  (/(r ).  An  example  of  this  good 
fit  IS  shown  in  fig.  6.  Note  that  the  assumption  of  F{x) 
being  a  single  exponential  decaying  function  implies 
that  P(  r(  /? ) )  is  of  Lorentzian  form  as  shown  in  eq. 
( 2  ).  This  function  has  a  diverging  second  moment, 
to  which  the  energy  density  in  the  fluid  is  propor¬ 
tional.  Therefore  G(t)  cannot  have  this  form  for  large 
values  of  We  indeed  observed  departures  from 
this  Lorentzian  form  for  P{y(R))  with  very  large 
values  of  I'lR)  (corresponding  to  very  small  f  for 
O  ( M  )  [10.11].  However,  these  observations  will  not 
be  discussed  further  here.  Most  theories  of  turbu¬ 
lence  concentrate  on  the  scaling  behavior  of  the  mo¬ 
ments  of  L(/?),  rather  than  in  P(,V(R))  itself.  Quite 
often.  Pi  I’iR) )  is  assumed  to  be  of  Gaussian  form. 
P( !'(/?))  -  exp{  -  (  f'f /?)/«(/?)  I but  this  form  of 
/*(  1  ■(/?))  is  clearly  contrary  to  our  findings. 

How  can  one  understand  that  the  exponent  C  in¬ 
creases  from  0  to  «  1  /3  as  the  fluid  becomes  increas¬ 
ingly  turbulent?  A  fundamental  understanding  of  this 
result  is  lacking,  but  it  can  be  said  that  the  observa¬ 
tion  IS  consistent  with  the  notion  that  the  turbulent 
active  region  is  a  fractal  [  12].  Let  the  fractal  dimen¬ 
sion  of  the  turbulent  region  be  D.  Since  the  turbulent 
energy  is  confined  to  active  regions  of  dimension 
D  <  3.  the  concentration  of  the  turbulent  energy  is  in¬ 
creased  to  smaller  regions,  relative  to  the  case  of  vol¬ 
ume-filling  turbulence.  Modifying  the  Kolmogorov 


Fig  6  A  t>pical  autocorrelation  Tunction  g{t)  versus  t  in  gnd 
How  The  solid  line  is  a  fit  to  the  incomplete  gamma  function. 


theory  to  take  this  effect  into  account  [13],  one  has 
u(R)~R'.  with  C=  ]( 1  +  Z)-3).  According  to  this 
model,  the  increase  of  C  from  0  to  *  i  corresponds  to 
an  increase  of  D  from  2  to  3. 

It  should  be  stressed  that  our  measurements  of  ^(  / ) 
described  above,  do  not  directly  give  information 
about  the  fractal  dimension  of  the  energy-containing 
eddies;  it  can  only  be  said  that  the  data  invite  such  an 
interpretation.  The  above  interpretation  of  the  data 
in  fig.  4  is  supported  by  the  recent  work  of  Shreeni- 
vasan  et  al.  [  3 ) .  They  measured  the  fractal  dimen¬ 
sion  of  the  interface  of  two  counter-flowing  fluids,  one 
of  which  has  been  dyed.  Such  measurements,  made 
in  the  vicinity  of  Rec,  support  the  conclusion  that  in¬ 
creasing  Re  above  Rec  increases  the  dimensionality 
of  the  turbulent  active  region.  With  one  adjustable 
parameter.  Rec.  the  data  of  Shreenivasan  et  al.  can  be 
directly  superimposed  on  the  measurements  in  fig.  4 
(3). 

Even  if  the  energy-containing  eddies  in  a  turbulent 
field  occupy  regions  with  dimensionality  less  than  3. 
it  is  not  necessary  that  the  entire  turbulent  region  be 
characterized  by  a  homogeneous  fractal.  Benzi  et  al. 
[14]  have  proposed  a  model  that  the  turbulent  re¬ 
gion  is  a  multifractal  object.  In  their  model  there  is  a 
probability,  x.  that  the  turbulent  region  is  space  fill¬ 
ing  (D=3)  and  a  probability.  1  -x.  that  D=2.  Our 
measurements  are  consistent  with  this  model,  pro¬ 
vided  one  makes  an  additional  assumption  that  x  is 
a  function  of  the  Reynolds  number.  The  details  of  the 
model  have  been  worked  out  for  a  general  function 
of  .rlRc)  and  fitted  to  experiment  [9].  At  present, 
however,  our  measurements  are  not  precise  enough 
to  confirm  that  a  multifractal  model  is  required  to 
explain  the  observations. 


4.  Concinding  remarks 

What  started  out  as  a  study  carried  out  in  the  time 
domain  ( the  measurement  of  f ) ).  has  ended  up  by 
yielding  spatial  information  about  turbulent  flow.  The 
homodync  experiments  provide  further  confirma¬ 
tion  of  the  notion  of  Mandelbrot  that  the  energy- 
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containing  eddies  are  fractal  in  their  geometrical 
structure.  This  finding  is  not  new.  What  seems  to  have 
gone  unnoticed  before,  is  that  the  fractal  dimension 
of  the  turbulence  changes  with  changing  Reynolds 
number,  when  some  critical  value  of  this  parameter 
is  exceeded.  The  interpretation  of  these  experiments 
makes  no  appeal  to  the  frozen  turbulence  assump¬ 
tion.  By  using  the  technique  of  photon  correlation 
homodyne  spectroscopy,  we  have  been  able  to  ob¬ 
serve  the  self-similar  behavior  of  turbulent  flows  at 
moderate  Reynolds  number  that  were  heretofore  re¬ 
garded  as  too  weak  to  exhibit  universal  features. 
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ABSTRACT.  Photon  Homodyne  Spectroscopy  (HCS)  and  Lziser  Doppler  Ve- 
locimetry  (LDV)  have  been  used  to  study  turbulent  velocity  fluctuations  V{R)  asso¬ 
ciated  with  eddies  of  size  R.  For  small  R  both  types  of  measurement  were  consistent 
with  a  model  in  which  the  active  regions  of  turbulence  lie  on  a  fractal  of  dimension 
D,  with  D  increasing  from  ~  2  to  ~  3  as  the  Reynolds  number  (Re)  increases  above 
some  threshold  value.  At  larger  eddy  sizes,  the  LDV  measurements  show  a  different 
scaling  of  the  velocity  fluctuations.  We  associate  these  larger  eddies  with  the  energy 
reservoir  that  feeds  the  inertial  cascade. 


INTRODUCTION 

It  is  known  that  the  motion  of  an  incompressible  fluid  is  regular  and  smooth  at  low 
Reynolds  number  (Re).  As  Re  increases,  the  flow  becomes  chaotic  and  irregular, 
and  finally  it  develops  into  tvirbxilence.  One  of  the  essential  aspects  of  turbulence 
is  that  the  large-scale  eddies  which  initially  contain  most  of  the  kinetic  energy,  spht 
into  smaller  s<^e  eddies,  and  the  process  is  repeated,  ag2dn  and  again,  until  the  size 
of  the  eddies  is  so  small  that  their  kinetic  energy  dissipates  into  heat.  Elementary 
theories  of  fully  developed  turbulence  are  based  on  the  hypothesis  that  the  small  scale 
statistics  of  turbvdent  flow  obeys  imiversal  scaling  properties,  such  as 

W'(R)  ~  R'’",  (1) 

where  V(R)  is  the  velocity  difference  associated  with  turbulent  eddies  of  size  R.  In 
general,  h  %  1/3  for  fully  developed  turbulence  (Kolmogorov  1941). 

We  describe  experiments  in  which  the  above  sw-similarity  is  observed  in  a 
weakly  turbulent  grid  flow,  even  though  the  maximum  Reynolds  numbers  reached  in 
these  experiments  is  not  very  large.  This  self-similarity  or  scaling  was  seen  in  two 
different  ways:  homodyne  photon  correlation  spectroscopy  (HCS)  and  laser  Doppler 
veiocimetry  (LDV). 


TECHNIQUES 

With  the  LDV  technique  (Durrani  zmd  Created  1977)  the  Doppler  shifted  light  scat¬ 
tered  by  small  particles  suspended  in  the  turbulent  fluid  is  mixed  with  the  reference 


beam.  The  detected  light  intensity  I{t)  is  therefore  modulated  by  the  Doppler  fre¬ 
quency,  q  •  v(r,  t),  where  q  is  scattering  vector  and  v(r,t)  is  the  local  velocity  of 
scatterers  in  the  fluid.  Our  interest  in  the  LDV  measurement  was  focused  on  the 
variance, 

<  V\r)  >=<  w(e  +  r)-  v(f')P  >  . .  (2) 

To  relate  <  V^(t)  >  to  the  spatial  velocity  variation,  <  V‘^{R)  >,  we  used  Taylor’s 
hypothesis  (frozen  turbulence  assumption)  which  assumes  R  =  C7t,  where  U  is  the 
mean  flow  velocity. 

With  the  HCS  technique  (Berne  and  Pecora  1976,  Tong  et  al.  1988),  one 
measures  the  intensity  correlation  function  g{t)  =<  I{t'  +  t)I(t')  >.  If  one  normalizes 
g{t)  so  that  g{t)  =  1-1-  const  •  G{t),  with  <j(0)  =  1  and  C?(oo)  =  0,  then 

G(f)  =  f^h{R)dR  r  dV{R)P{V{R),R]co3(qV{R)t).  (3) 

Jo  J —oo 

Here,  V{R)  is  the  component  of  the  instantaneous  velocity  difference  V(jR)  along  the 
scattering  vector  q,  and  P{V{R).  R)  is  its  probability  density  function.  The  scattering 
volume  viewed  by  a  photodetector  is  assumed  to  be  quasi-one-dimensional  with  length 
L,  and  h{R)  =  2(1  —  R/L)/L  is  the  probability  density  of  finding  a  pair  of  particles 
separated  by  a  distance  R  in  the  scattering  volume.  Elquation  (3)  means  that  the 
li(^t  scattered  by  each  pair  of  particles  contributes  a  phase  factor  cos{qV{R)t)  (due 
to  the  frequency  beating)  to  the  intensity  correlation  function  gU),  and  g{t)  is  the 
incoherent  sum  of  these  ensemble  averaged  phase  factors  over  all  the  particle  pairs  in 
the  scattering  volume. 


EXPERIMENTAL 

The  fluid  flow  was  generated  in  a  closed  water  tunnel  with  a  variable  speed  pump. 
Undesirable  flow  disturbances  produced  by  the  pump  and  piping  system  were  damp^ 
out  in  the  plenum  chamber  on  the  high  pressure  side  of  the  grid.  The  grid,  which 
generates  the  turbulence,  consists  of  an  array  of  rods  of  width  0.016  cm  and  with  a 
mesh  size,  M  =  0.85  cm.  The  Reynolds  number  is  defined  as  Re  =  U M/ u,  where  v 
is  the  kinematic  viscosity  of  water. 

The  tunnel  itself  (Pak  et  al.  1991)  is  110  cm  long  and  its  square  cross  section 
is  10  cm  X  10  cm.  Polished  plexiglass  windows  are  furnished  on  four  sides  of  the  water 
tunnel  to  admit  the  incident  laser  beam  and  to  observe  the  scattering.  A  recovery 
tank  is  located  at  the  low  pressure  end  of  the  tunnel.  Its  large  volume  slows  the  flow 
before  the  water  returns  to  the  pump  intake.  The  water  was  distilled  and  seeded 
with  polystyrene  spheres  of  uniform  size.  These  particles  scatter  the  light  for  both 
the  HCS  and  LDV  measurements.  The  observation  point  for  both  experiments  was 
on  the  aixis  of  the  tunnel  and  21  cm  down  stream  from  the  grid,  where  the  spanwise 
velocity  profile  is  fiat.  All  measurements  were  made  at  room  temperature. 

In  the  HCS  experiment  a  mildly  focused  argon-ion  laser  beam  entered  the  test 
section  from  above,  md  the  scattering  measurements  were  made  in  the  horizontal 
plane.  A  lens  was  placed  in  such  a  position  as  to  form  an  1:1  image  in  the  plane  of 
an  adjustable  slit.  It  is  the  slit  width  L  which  controls  the  scattering  volume  size 
viewed  by  the  photomultiplier.  For  the  LDV  measxirement  we  used  standard  tech¬ 
niques  for  measuring  the  stream-wise  component  of  the  velocity  v^{t).  The  Taylor’s 
hypothesis  is  expected  to  be  valid  only  if  the  rms  velocity  fiuctuations  aroimd  the 
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Figure  1.  Measured  G(t}'s  for  (a)  the  turbulent  grid  flow  of  Re  =  1800  and  L  = 
0.7  mm  and  (b)  turbulent  Couette  flow  of  Re  =  8.3  x  10*  and  L  =  1.0  mm. 


mean,  (vg(t)  —  C/)^  >,  are  small  compared  to  U  itself.  For  all  the  measurements 
reported  here,  (v*(f)  —  17)~  >/U  was  less  than  4%. 


RESULTS 

In  the  HCS  measurement,  it  was  found  that  over  a  limited  range  of  slit  widths  and 
Reynolds  numbers,  G{t)  has  the  scaling  form,  G(t)  oc  G(gu(L)t).  Figure  1(a)  shows 
a  typical  G(t)  in  the  water  tunnel.  Here  Re  is  1800  and  the  slit  width  is  0.7mm.  We 
define  u(L)  from  the  inverse  of  the  decay  time  of  G(t),  i.e.  l/T  =  qu(L).  Therefore 
u(L)  is  the  characteristic  turbulence  velocity  difference  at  length  sc^e  L.  The  above 
scaling  form  of  G(t)  was  seen  only  when  Re  exceeded  some  critical  value,  Re^  ~  300. 
It  was  also  found  that  u(L)  has  a  scaling  form,  u{L)  ^  L°.  As  Re  increases  from 
Re  :::  300  to  800,  a  increases  from  0  to  1/3  (the  Kolmoeorov  value),  and  beyond 
Re  ~  800.  a  saturates  at  the  value  of  1/3  (see  Fig.  2(a)).  The  minimum  usable  value 
of  L  is  limited  by  the  laser  beam  width  to  0.2  mm.  The  upper  limit  of  L  was  set  by 
optical  coherence  effects  (Berne  and  Pecora  1976)  to  approximately  0.2cm. 

Figure  3  is  a  log-log  plot  showing  a  typical  <  V^{R)  >  vs  R  of  LDV,  sub¬ 
tracting  R- independent  contributions  from  the  instrumental  phase  fluctuations  (Pak 
et  al.  1991).  These  measurements  were  made  at  Re  =  330  (curve  (a))  and  Re  =  1400 
(curve  (Iri)  respectively.  The  curves  have  been  shifted  vertically  to  improve  visibility. 
In  the  LDV  measurements  the  smallest  R  was  limited  by  the  scattering  volume  size 
0.3  mm,  and  in  the  opposite  limit  of  R  is  approximately  the  mesh  size,  M  =  0.85  cm. 
For  Re  ^  300(Rce),  there  we  two  scaling  regimes  in  which  <  V^{R)  >~  R^  for 
R^O.lScm  and  <  V^(R)  >~  RP  for  R^  0.15cm.  Curves  (b)  and  (c)  of  Fig.  2 
are  scaling  exponents  C/2  for  small  R  and  77/2  for  large  R  respectively.  For  small  R, 


Figure  2.  Scaling  exponents  in  the  inertial  subrsoige  and  the  energetic  subrange 
measured  by  LDV  and  by  the  HCS  techniques:  (a)  a{Re)  in  the  inertial  subrange 
measured  by  HCS.  (b)  (/2(Re)  in  the  inertial  subrange  measured  by  LDV.  (c)  T)/2{Re) 
in  the  energetic  subrange  measured  by  LDV. 


where  both  techniques  probe  overlapping  ranges  of  R,  as  a  ingreftsfis  from  0  to  1/3.  C/2 
decreases  from  1/2  to  1/3.  Both  exponents,  ct  and  C/2  saturate  at  1/3  for  Rs^  800. 
For  large  R,  as  Re  increases,  tj/2  decreases  from  0.8  and  saturates  at  0.5.  Below  we 
offer  a  possible  explanation  as  to  whv  q  and  C  should  vary  oppositely  with  increasing 
Re. 


DISCUSSION 

To  interpret  the  LDV  and  HCS  measurements  at  small  R,  we  invoke  the  (3  model 
(Mandelbrot  1976,  Frisch  et  al.  1978),  According  to  this  model  the  turbident  kinetic 
energy  in  the  inertial  range  lies  in  an  “active”  region  having  fractal  dimension  D. 
With  this  model  P{V{R),R)  can  be  factorizable  to  the  form.  P*(tt(i?),i2)  •  i?“^. 
Here  P*{u{R),R)  is  the  probability  density  of  u(H)  in  the  active  region  and  the 
factor  i2®“®  is  proportional  to  the  sparcity  of  active  turbulent  regions  at  small  R. 
Therefore  the  product  of  R^~^  and  h{  R)  is  proportional  to  the  probability  of  finding 
a  particle  pair  in  an  active  region  within  the  slit.  These  factors  have  little  effect 
on  u{R).  For  this  reason,  an  HCS  experiment  of  u(I)  is  a  measurement  of  velocity 
differences  in  the  active  regions  of  the  turbulence  only.  Prom  this  model  one  readily 
derives  that  a  =  (Z)  —  2)/3.  Therefore  an  increase  of  a  from  0  to  1/3  implies  that  D 
is  increasing  from  2  to  3. 

The  LDV  technique  for  small  R,  unlike  HCS,  probes  both  the  “active”  and 
“inactive”  regions  of  the  flow,  since  both  regions  are  moved  past  the  observation 
point  at  the  mean  flow  speed  U.  Because  of  the  sparcity  of  active  turbulent  regions, 
<  \^{R)  u^R)  R^^  ~  R^  with  C  =  1-ct  =  (5-D)/3.  Accordingly,  C/2  should 


R  (cm) 

Figvire  3.  Log-log  plot  of  <  V^{R)  >  vs  R  from  the  LDV  measurements.  Shown 
here  are  measurements  at  Re  =  1400  (a)  and  Re  =  330  (b).  The  daished  and  solid 
lines,  of  slope  ^/2  and  77/2  respectively,  are  discussed  in  the  text. 


decrease  from  1/2  to  1/3  as  a  increases  from  0  to  1/3.  The  small* i2  range  (inertial 
range),  where  the  above  fractal  model  can  be  applied,  is  limited  from  bdow  by  the 
Kolmogorov  dissipation  scale.  The  inertial-range  eddies  are  fed  by  a  reservoir  of 
larger,  energy-containing  eddies  which  normally  are  not  associated  with  self- similarity. 
Nevertheless  the  measurements  of  <  V^{R)  >  at  large  R  (the  soHd  Unes  of  Fig.  2 
and  curve  (c)  in  Fig.  3)  seem  to  define  an  exponent  (77)  associated  with  these  energy- 
containing  eddies.  Since  the  range  of  R  where  77  is  well  defined  is  much  less  than  a 
decade,  ri{Re)  may  have  little  physical  meaning. 

Recently  we  have  commenced  HCS  measurements  in  a  flow  between  two  con¬ 
centric  cylinders  (Couette  Cell),  with  the  inner  one  rotating  (Tong  et  al.  1990).  When 
this  flow  was  highly  turbulent,  G(t)  (and  by  implication  of  P{V{R),  R))  was  clearly 
of  different  functional  form  than  for  the  water  tunnel  measurements.  Curve  (b)  of 
Fig.l  is  G{t)  for  the  Couette  flow  at  Re  =  8.3  x  10*  and  slit  width  L  =  Ittitti  when 
the  laser  beam  passed  through  the  gap  radially. 


CONCLUSION 

The  combined  LDV  and  HCS  experiments  reported  here  provide  additional  evidence 
that  in  the  inertial  range  of  eddy  sizes,  the  fractal  dimension  D  of  the  turbulence 
increases  with  increasing  Reynolds  number.  The  observations  are  consistent  with  an 
increase  in  D  from  ~  2  to  ~  3  as  Re  increases  above  some  system-dependent  threshold 
value.  Rce.  This  assertion  is  based  on  the  recognition  that  these  two  measurement 
schemes  average  over  the  velocity  fluctuations  V{R)  in  different  ways.  The  HCS 
technique  probes  the  active  region  only,  whereas  the  LDV  meastirements  average 
over  the  inactive  regions  as  well. 
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Homodyne  photon  correlation  spectroscopy  (HCS)  and  laser  Doppler  velocimetry 
(LDV)  were  used  to  study  the  probability  density  function  of  velocity  differences 
6v(l)  between  points  in  the  fluid  separated  by  a  distance  t.  These  two  measuring 
schemes  yield  different  results  for  the  probability  density,  P{6{v{t))dSv{i).  Both 
types  of  measurement  show  that  P(Sv(f))  is  well  approximated  by  the  product  of 
Gaussian  and  Lorentzian  factors,  with  P  being  more  Lorentzian-like  for  the  HCS 
measurements  and  more  Gaussian-like  for  the  LDV  measurements.  The  HCS  method 
probes  the  spatial  fluctuations  directly,  whereas  the  LDV  technique  records  temporal 
fluctuations  in  velocity  and  relates  them  to  the  spatial  fluctuations  through  the  Taylor 
hypothesis. 

PACS  numbers:  47.25.Ae,  47.25. Cg 
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In  turbulent  flow  the  quantity  of  interest  is  not  so  much  the  velocity  itself,  v(r,  t), 
as  the  instantaneous  difference  in  velocity  at  two  points  separated  by  a  distance  £. 
This  quantity  Sv(£)  =  t;(r  +  £,f)  —  v(r,t)  is  characterized  by  a  probability  density 
function,  P(Sv(£))  dSv(£),  where  t;  is  the  measured  component  of  v.  Herein  we  report 
measurements  of  P(Sv(£))  in  turbulent  flow  generated  by  a  grid  in  a  water  tunnel.  The 
experiments  explore  a  range  of  ^-values  and  span  a  range  of  Reynolds  numbers  under 
conditions  that  possibly  lie  between  chaotic  behavior  and  fully  developed  turbulence. 

Our  measurements  of  P(6t;(£))  will  be  compared  with  a  measured  probability 
density  function  P'(6u(f))  deduced  from  a  record  of  the  velocity  as  a  function  of  time, 
recorded  at  a  single  point  in  the  turbulent  fluid.  Both  types  of  measurement  were 
made  at  almost  exactly  the  same  position  in  the  turbulent  stream.  To  deduce  spatial 
information  from  the  time  record  of  the  the  velocity,  it  is  necessary  to  invoke  the 
so-called  Taylor  hypothesis  (frozen  turbulence  assumption)  [1],  which  simply  means 
that  v(<)  is  replaced  by  v{xfU)  where  U  is  the  mean  flow  velocity.  Virtually  all 
measurements  of  the  probability  density  function  for  velocity  differences  on  scales  £ 
have  been  made  in  this  way  [2].  It  is  this  probability  density  which  we  have  designated 
as  F. 

In  the  work  reported  here,  P  and  P'  are  found  to  be  very  different.  As  for  P,  it  is 
well  represented  by  the  product  of  Gaussian  and  Lorentzian  terms.  The  Lorentzian 
factor  is  found  to  have  a  much  smaller  width,  U[,{£),  than  the  width  UG(£)y  of  the 
Gaussian  factor.  The  new  measurements  reported  substantiate  those  made  in  a  much 
smaller  water  tunnel  [3-5].  The  function  P*  is  much  more  Gaussian  in  form  than  P, 
but  it  too  gives  a  heayier  weighting  to  large  velocity  fluctuations  than  does  a  purely 
Gaussian  function  or  even  an  exponential  one.  If  one  parameterizes  P'  by  the  product 
of  Gaussian  and  Lorentzian  factors,  the  widths  of  these  two  factors  are  comparable 
in  magnitude.  For  our  measurement  at  least,  the  frozen  turbulence  assumption  fails. 
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i.e,  P  and  P*  are  conveying  different  information  [6]. 

The  schemes  for  measuring  P'  and  P  are  laser  Doppler  velocimetry  (LDV)  and 
homodyne  photon  correlation  spectroscopy  (HCS),  respectively.  The  former  technique 
is  a  standard  one  [7],  whereas  the  HCS  method  has  been  used  only  recently  [3,8-11]. 
Both  of  these  methods  require  the  seeding  of  the  fluid  with  small  particles  which  follow 
the  local  flow  and  scatter  light.  The  HCS  method  is  sensitive  to  velocity  differences 
rather  than  the  velocity  itself,  because  it  registers  the  beating  of  Doppler  shifts  coming 
from  pairs  of  particles  moving  relative  to  each  other  [3].  As  will  be  discussed  below, 
the  optical  scheme  used  here  records  velocity  differences  for  all  values  of  £  from  the 
smallest  eddy  size  present  out  to  a  size  L,  which  is  the  width  of  a  slit  through  which 
the  light  passes  before  it  reaches  the  photodetector.  Because  one  measures  a  scattered 
intensity  rather  than  a  scattered  electric  field,  the  HCS  method  yields  information 
about  the  symmetric  part  of  the  probability  density  only.  It  is  this  part  of  the  full 
probability  density  function  that  we  label  as  P{6v{£)).  Past  measurements  of  P'  in 
highly  turbulent  flows  show  that  this  function  is  not  always  symmetric  [12, 13].  In 
the  present  LDV  measurements  P'{6v{())  is  also  asymmetric,  though  the  asymmetry 
is  small. 


Our  HCS  measurements  of  P{Sv{£))  show  it  to  be  well  approximated  by  the 
product  of  Lorentzian  and  Gaussian  functions  only  if  the  Reynolds  number  of  the 
flow.  Re,  exceeds  some  critical  value,  Reg  ^  300.  Here  we  define  Pc  as  Pe  =  aU/v, 
where  a  is  the  mesh  size  of  the  grid  (a  =  8.5  mm)  and  i/  is  the  kinematic  viscosity 
of  water  (O.Ol  cm^/sec).  The  widths  of  the  Lorentzian  and  Gaussian  factors  are 
explicitly  defined  by  the  equation. 


ezp(-uc(£)^/2uG(£n 

rer/c(ui(£)/y/2uG(£)) 


exp(6v{£f  I2ug{£?) 


UL(f) 

Sv(£y-i-uc(£y' 


(1) 


The  ratio  of  these  widths.  A/  =  uc(f)/u£(f)  is  found  to  be  approximately  4.  This 
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ratio  is  only  weakly  dependent  on  i  and  Rt  in  the  range  of  parameters  probed  in  this 
work  namely,  0.3  mm  <  Z  <  2.5  mm,  300  <  Re  <  1800. 

As  discussed  below,  the  quantity  actually  meaisured  with  the  HCS  scheme  is  the 
intensity  correlation  function,  of  light  scattered  by  the  seed  particles.  This 
function  is  observed  to  have  a  self-similar  form,  which  it  would  not  have  unless  both 
widths,  and  ug(^)  are  proportional  to  with  a  common  value  of  the  exponent 
h.  In  this  experiment,  as  in  the  earlier  work,  using  a  smaller  water  tunnel,  h  is 
observed  to  increase  with  increasing  Re  from  h  ~  0  to  h  ~  1/3,  the  Kolmogorov 
value  [14]. 

The  experiments  were  carried  out  in  the  water  tunnel  that  was  10  cm  x  10  cm  in 
cross  section  and  1.1  m  in  length.  The  measurements  were  made  at  an  axial  point  25 
cm  downstream  from  the  grid.  The  spanwise  prohle  of  the  flow  was  flat,  making  the 
time-averaged  velocity  gradient  at  the  measuring  point  for  P  and  P'  ignorably  small 
(11). 

The  water  was  seeded  with  polystyrene  particles  (diameters  0.106 that  were 
small  enough  to  follow  the  local  flow.  The  light  source  was  the  mildly  focused  beam 
from  an  argon  ion  laser  operating  the  wavelength  A  =  488  nm.  The  incident  beam 
travelled  perpendicular  to  the  flow  direction.  The  scattered  light  was  detected  in  a 
direction  that  was  also  perpendicular  to  the  flow.  The  incident  beam  traced  out  a 
clearly  visible  thin  line  in  the  flowing  water,  the  diameter  of  this  beam  being  less 
than  0.1  mm.  The  scattered  light,  at  a  scattering  angle  9  s  90**,  was  imaged  at  1:1 
magnification  on  a  slit  of  adjustable  width  Z.  This  width  is  a  crucial  parameter  in 
the  experiment  because  it  determines  the  maximum  measurable  eddy  size,  (  —  L, 
to  which  g{t)  is  sensitive.  More  than  1.5m  behind  the  slit  was  a  photomultiplier 
which  recorded  the  scattered  intensity,  I{t).  The  output  of  the  photomultiplier  is  a 
train  of  identically  shaped  pulses  which  are  sent  to  the  correlator,  whose  output  is 
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the  intensity  correlation  function  g{t)  =<  I{t*  +  >  /  <  I{t')  >^. 

Under  approximations  which  are  well  satisfied  in  these  experiments,  one  can  show 
[3]  that  g{t)  is  given  by 

g{t)  =  1  +  Ait)  £  b{i)dl P  Pi6v{t))cosik6vit)t)d8v{t))  =  1  -h  G(t),  (2) 

where  the  scattering  vector  k  is  of  magnitude  k  =  2nkoSin{6f2),  and  6vi()  is  the 
component  of  6v{£)  along  the  direction  of  k.  The  refractive  index  of  the  scattering 
medium  (water)  is  n  =  1.33.  Quite  generally,  G(/)  is  a  decaying  function,  its  char¬ 
acteristic  decay  time  that  is  of  the  order  of  l/jku(L),  where  u{L)  is  the  characteristic 
velocity  difference  over  a  slit  width  L.  The  function  6(£)  =  (2/Z,)(l  —£/!,)  is  the  prob¬ 
ability  that  a  pair  of  particles,  separated  by  £,  are  to  be  found  in  the  sht  [3, 11).  This 
form  of  b{()  is  correct  only  if  L  is  much  greater  than  the  beam  diameter  [3,8].  The 
factor  A{t)  represents  the  Brownian  diffusion  of  the  seed  particles  and  takes  account 
of  the  fact  that  even  in  the  absence  of  turbulent  flow,  G{t)  will  decay  slowly  and  is  of 
the  form  Ait)  —  exp[-2Dk%  where  D  is  the  diffusion  constant.  This  parameter  is 
determined  from  the  diameter  of  the  seed  particle  and  the  viscosity  of  water  (15).  All 
of  the  experimental  data  have  been  corrected  for  this  Brownian  motion  factor  [llj. 

According  to  Eq.  (2),  if  P(^u(£))  has  the  simple  scaling  form  P(6u(£))  = 
QiSvi())/uii)),  and  if  ui()  ~  then  Git)  has  the  scaling  form,  Gik,t,L)  =  G(x), 
where  x  =  kuii)t.  Therefore  a  plot  of  log  Git)  vs  log  t,  for  fixed  Re  but  various  values 
of  L,  should  produce  curves,  all  of  which  should  be  superimposable  when  translated 
along  the  time  axis.  This  is  indeed  what  was  observed  in  the  earlier  experiments  [3] 
and  in  the  present  ones.  Figure  1  shows  G(t)  measured  at  the  three  slit  widths  L  = 
0.5  mm,  1.1  mm,  and  2  mm,  the  Reynolds  number  being  Re  =  1490.  The  function 
G(0  has  been  normalized  so  that  it  approaches  unity  in  the  limit  t  — »  0  and  zero  as 
<  — ►  00.  Having  confirmed  that  P  is  self-similar,  we  now  examine  its  functional  form. 
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Figure  2  is  a  semi-log  plot  of  G(t)  at  Jie=  1490,  Zr  =  1.3  mm.  Also 
shown  is  the  best  fit  to  Eq.  (2),  assuming  that  P  is  purely  Gaussian  (P  = 
ll/\/^uo(^)]exp(—Sv(£y/2uG(£y))  (solid  line),  and  purely  Lorentzian  (P  = 
^l(£)/^I^l(^)^  +  (broken  line).  These  two  lines  correspond  to  ul(L)  = 

ua(L)  =  0.216  cm/sec.  Variation  of  these  parameter  could  not  produce  a  good  fit 
to  these  data  with  either  the  Gaussian  or  Lorentzian  probability  density  function, 
though  the  Lorentzian  fit  is  clearly  superior.  An  exponential  form  for  P  is  even  more 
at  odds  with  our  observations,  though  it  is  observed  in  other  experiments  [12].  We 
have  shown  a  Gaussian  form  for  P(^i;(£))  because  it  is  sometimes  observed  for  large 
values  of  £  [16].  A  probability  density  of  dominantly  Lorentzian  form  gave  a  satisfac¬ 
tory  fit  to  the  measurements  in  Ref.  [4].  We  defer  a  discussion  of  other  suggested  and 
observed  forms  of  P(6u(^)). 

To  proceed  further  we  take  note  of  the  fact  that  the  cosine  factor  in  Eq.  (2)  implies 
that  G(t)  must  approach  the  origin  with  zero  slope  (G(t)  l  —  const.P),  even  though 
experimental  limitations  on  the  electronics  bar  an  unambiguous  observation  of  this 
behavior.  We  therefore  attempted  the  fit  to  G(t)  with  the  product  of  a  Lorentzian 
and  a  Gaussian  function,  the  widths  of  these  factors  being  adjustable  parameters. 
As  seen  in  the  inset  of  Fig.  2  (solid  line),  it  was  possible  to  obtain  a  very  good  fit 
to  the  data  with  this  form  of  P.  A  similarly  good  fit  with  this  form  of  G(t)  was 
obtained  at  all  values  of  L  and  Be.  For  each  of  the  parameter  pairs,  ug(^)  and  ui(£), 
we  have  evaluated  M{£,Re)  and  find  it  to  be  approximately  4,  increasing  slightly 
with  decreasing  slit  width  (At  Re  =  1490,  A/  =  2.7  at  Z  =  1.3  mm,  and  M  —  3.6 
at  Z  =  0.5  mm).  There  is  no  discernable  trend  with  changing  Re.  Similar  values 
of  M  were  obtained  in  Ref.  [4].  Onuki  [17]  has  suggested  that  these  two  widths 
characterize  velocity  fluctuations  inside  and  outside  of  the  active  regions.  Abundant 
experiments  suggest  that  the  active  regions  of  the  turbulent  dissipation  lie  on  a  fractal 
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or  multi-fractal  (18]. 

It  is  instructive  to  compare  our  measurements  of  G(t)  with  a  prediction  that  P 
is  a  log-normal  distribution,  since  this  functional  form  of  P  has  a  credible  theoretical 
basis  and  many  experiments  have  been  aimed  at  testing  the  log-normal  theory  [2, 
16, 19).  Here  we  compare  our  mcMurements  with  a  modified  log-normal  theory  due 
to  Castaing  et  al.  [13].  This  theory,  which  is  free  of  iin  inconsistency  contained  in 
the  original  log-normal  theory  [2],  successfully  explains  recent  measurements  of  the 
temperature- difference  density  function  P(6T{£))  as  measured  in  Rayleigh  Benard 
convection  by  the  Chicago  group  [20].  Castaing  assumes  that  6v{i)  has  a  Gaussian 
distribution  with  variance  u(£),  but  that  u{£)  is  a  random  variable  with  a  log-normsd 
distribution.  The  specific  form  of  P  is 

Px{6v{£))  =  f  exp{—6v{£)^l2u^)ezp{—lTp{u/u{£))/2X^)dulv},  (3) 

Jo 

where  u(£)  is  the  most  probable  variance  of  6v{£)  and  A  is  the  variance  of  In  u.  To 
compare  this  prediction  with  our  measurements,  we  assume  that  u(£)  may  be  taken  to 
be  the  measured  value  of  ul{L)  (since  the  Lorentzian  f2u;tor  dominates  in  G{t)).  The 
parameter.  A,  is  taken  to  have  the  form  proposed  by  Castaing,  namely  A  =  {£/£o)^, 
with  0  =  0.24.  The  parameter  £„  is  set  equal  to  the  grid  size  of  8.5  mm.  Inserting 
the  above  expression  for  P  into  Eq.  (2)  yields  a  G{t)  which  is  plotted  as  the  broken 
line  in  the  inset  of  Fig.  2.  We  note  that  a  log  normal  form  of  P,  such  as  in  Eq.  (3), 
assures  that  G(<)  will  not  have  the  self-similar  form  so  clearly  seen  in  Fig.  1.  With 
the  above  values  of  and  u(£),  the  log-normal  theory  is  clearly  not  a  good  fit  to  the 
data  in  Fig.  2. 

Next  we  turn  to  the  LDV  measurements  of  P'(5u(£)).  With  the  LDV  technique, 
the  seed  particles  traverse  optical  interference  fringes,  modulating  the  intensity  of  the 
scattered  light  received  by  the  photomultiplier.  A  commercially  available  counter- 
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signal  conditioner  records  the  frequency  of  these  pulses  and  hence  determines  v{t)  at 
the  observation  point.  The  Taylor  hypotheses  is  then  invoked  to  convert  the  v{t)  into 
a  temporally  fluctuating  signal,  t;(x),  where  z  is  a  coordinate  in  the  flow  direction  [7]. 
The  Taylor  hypothesis  is  generally  assumed  to  be  valid  when  the  ratio  /  of  the  rms 
fluctuations  in  the  velocity  about  its  mean  value,  U,  is  a  small  fraction  of  U  itself. 
This  condition  was  well  satisfied  in  the  present  experiments  where  /  was  always  less 
than  0.04.  By  splitting  up  the  time  record,  v(t),  into  segments  of  equal  length  and 
replacing  6i  by  tjU,  we  create  an  the  ensemble  of  velocity  differences,  ^v(^),  from 
which  P'{6v{i))  d6v{6C)  can  be  constructed  with  a  histogram. 

Fig.  3  is  a  semi-log  plot  of  an  LDV  measurement  of  P'  vs  at  i2c  =  1490 

and  ^  =  1.3  mm.  The  solid  line  in  the  figure  is  a  best  fit  to  these  data,  assuming  that 
P*  is  of  the  Lorentzian-Gaussian  form  of  Eq.  (1)  (If  P'  were  a  Gaussian  function,  the 
data  points  would  lie  on  a  straight  line  in  this  figure).  The  parameters  producing  this 
fit  are  ul(L)  =  0.32  cm/sec  and  ua{L)  =  0.47  cm/sec,  giving  M  =  1.4.  This  ratio  is 
much  smaller  than  M  obtained  from  the  HCS  experiments.  Those  LDV  measurements 
which  spanned  the  same  range  of  Re  and  L  as  the  HCS  experiments,  could  be  fitted 
with  a  density  function  containing  the  product  of  Lorentzian  and  Gaussian  factors. 
However,  the  exponent  h{Re)  could  not  be  determined  without  making  questionable 
assumptions  [11]. 

The  broken  line  in  Fig.  3  is  the  function  P  extracted  from  the  HCS  measurements 
made  using  the  same  values  of  Re  and  L  as  the  LDV  data  in  the  figure.  The  param¬ 
eters  associated  with  this  HCS  measurement  are  u^(L)  =  0.26  cm/sec  and  ua{L)  = 
0.69  cm/sec.  Note  that  the  Gaussian  factor  is  weighted  much  more  heavily  than  for 
the  HCS  measurements,  i.e.  M  is  smaller.  The  very  different  shape  of  the  two  curves 
in  Fig.  3  leaves  no  doubt  as  to  the  failure  of  the  frozen  turbulence  assumption  here. 
Both  P  and  P'  give  a  greater  weight  to  velocity  fluctuations  of  large  magnitude  than 
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either  a  Gaussian  or  exponential  function. 

The  failure  of  the  frozen  turbulence  assumption  indicates  that  the  small-eddies, 
which  are  being  probed  here,  do  not  maintain  their  geometrical  shape  as  they  are 
convected  past  the  observation  point  at  the  mean  velocity,  U.  A  possible  cause  is 
the  coupling  between  velocity  fluctuations  of  many  different  sizes,  a  point  that  has 
received  much  theoretical  attention  [21,22]. 

In  summary,  we  have  carried  out  measurements  which  expose  the  Taylor  hypoth¬ 
esis  (frozen  turbulence  assumption)  to  an  experimental  test  and  have  measured  the 
functional  form  of  the  probability  density  P(Sv{£))  over  a  range  of  eddy  sizes  £  in 
the  Reynolds  number  range  300  to  1800.  The  function  P  is  well  represented  by  the 
product  of  Gaussian  and  Lorentzian  factors  but  not  by  either  factor  alone  or  by  ex¬ 
ponential  or  log-normal  functions.  These  observations  indicate  the  importance  of 
probing  turbulence  in  both  the  time  and  space  domains. 
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FIGURES 


FIG.  1.  Superimposed  log-log  plot  of  G(t)  for  the  three  indicated  slit  widths,  £  and 
at  the  Reynolds  number  indicated. 


FIG.  2.  Semi-log  plot  of  G(i)  at  Re  =  1490,  L  =  1.3  mm.  The  solid  line  is  a  Lorentzian 
lit  to  the  data,  and  the  broken  line  is  a  Gaussian  fit.  The  characteristic  widths  vi{L)  and 
uc;(£)  of  both  factors  is  0.216  cm/sec.  The  inset  shows  the  same  data  fitted  to  £q.  (1) 
(solid  line)  and  to  the  log-normal  model  (broken  line)  discussed  in  the  text. 


FIG.  3.  Semi-log  plot  of  the  probability  density  function  as  measured  by  laser  Doppler 
velocimetry,  the  abscissa  being  6v{l)^.  The  solid  line  is  a  fit  to  these  data  using  Eq.  (1)  and 
the  broken  line  is  the  function  P  extracted  from  the  HCS  measurements,  using  Eqs.  (1)  and 
(2).  The  measurements  were  at  ResHOO  and  t  sl.3  mm. 
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